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Abstract 

The spin fractions and deep inelastic lepton structure functions of the proton 
are analyzed using chiral field theory involving Goldstone bosons. A detailed 
comparison with recent chiral models sheds light on their successful descrip- 
tion of the spin fractions of the proton as being due to neglecting helicity 
non-flip chiral transitions. This approximation is valid for zero mass but not 
for constituent quarks. Since the chiral spin fraction models with the pure 
spin-flip approximation reproduce the measured spin fractions of the proton, 
axialvector constituent-quark-Goldstone boson couplings are found to be in- 
consistent with the proton spin data. Initial quark valence distributions are 
then constructed using quark counting constraints at Bjorken x — >1 and Regge 
behavior at x — > 0. Sea quark distributions predicted by chiral field theory 
on this basis have correct order of magnitude and shape. The spin fractions 
also agree with the data. 
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I. INTRODUCTION 



The nonrelativistic quark model (NQM) explains qualitatively many of the strong, elec- 
tromagnetic and weak properties of the nucleon and other octet (and decuplet) baryons in 
terms of three valence quarks whose dynamics is motivated by quantum chromodynamics 
(QCD), the gauge field theory of the strong interaction. The effective degrees of freedom 
at the scale Aqcd are dressed or dynamical quarks along with Goldstone bosons which are 
expected to emerge in the spontaneous chiral symmetry breakdown (xSB) of QCD that is 
characterized by quark condensates. While the NQM description of baryon states has not 
been derived from first principles of QCD, solutions of Schwinger-Dyson equations for light 
quarks with models for color confinement lead to a momentum dependent quark mass m q (p 2 ). 
Such dynamical quarks become constituent quarks when m q (p 2 ) is approximated by a con- 
stant m q (0) at the scale Aqcd- In Nambu-Jona-Lasinio models |]J of the x^B patterned 
after superconductivity and driven by a zero-range scalar and pseudoscalar quark-quark 
interaction, the cut-off A^jl ~ 631 MeV of the effective chiral field theory provides an esti- 
mate for the momentum scale where the dynamical quark mass from the gap equation has 
fallen to the current quark mass. A more realistic treatment of color confinement is provided 
by the dual superconductor picture in the 't Hooft abelian gauge @ where QCD monopoles 
generated by the abelian gauge fixing condense in the QCD vacuum confining quarks by 
the dual Meissner effect. Here the dynamical xSB leads to a quark mass consistent with 
the QCD running whose monopole shape parameter A m ~ Aq C dVW ~ 630 MeV is close to 
A njl ■ At this momentum scale current quarks become the relevant degrees of freedom. 

At momentum scale Aqcd chiral perturbation theory (xPT allows incorporating 
systematically the chiral dynamics of QCD. Chiral field theory applies to scales from Aqcd 
up to a (presumed) chiral symmetry restoration scale A x = 4iif n = 1169 MeV where /„- = 93 
MeV is the pion decay constant. In this scenario, inside hadrons at distances that are 
smaller than the confinement scale Aq C d, but larger than distances where perturbative 
QCD (pQCD) applies, quarks and Goldstone bosons are the effective degrees of freedom. 
In principle, chiral field theory dissolves a dynamical or constituent quark into a current 
quark and a cloud of virtual Goldstone bosons, but here only one layer of Goldstone bosons 
is included. Thus, the effective degrees of freedom and interactions of chiral field theory are 
consistent with those of chiral perturbation theory. 

Chiral fluctuations of valence quarks inside hadrons, q^i — > qi,-\ + (qq')o, into pseudoscalar 
mesons, (qq')o, of the SU(3) flavor octet of Goldstone bosons, were first applied to the spin 
problem of the proton in ref . || . It was shown that chiral dynamics can help one understand 
not only the reduction of the proton spin carried by the valence quarks from AS = 1 in 
the NQM to the experimental value of about 1/3, but also the reduction of the axial vector 
coupling constant from the NQM value 5/3 to about 5/4. It is well known that relativistic 
effects reduce the axial charge further, and this causes problems for the spin fractions in 
relativistic quark models || forcing the proton size parameter to a rather large value of 
about 1 fm. This problem will no longer occur when we replace the NQM by constraints 
from the SU(6) spin-flavor wave function of the proton. In addition, the violation of the 
Gottfried sum rule GSR = 1/3 + 2{u — d)/3, which signals an isospin asymmetric quark 
sea in the proton, became plausible. The SU(3) symmetric chiral quark model explains 
several spin and sea quark observables of the proton, but not all of them. The data |8]|| call 
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for SU(3) breaking because some of the spin fractions such as the ratio of triplet to octet 
spin fractions A 3 /A 8 =5/3 and the weak axial vector coupling constant of the nucleon, 
9^a = F + T>, disagree with experiments in the SU(3) symmetric case. In [|T0| , [TT]1 the effects 
of SU(3) breaking were more systematically built into these chiral models and shown to 
lead to a remarkable further improvement of the spin and quark sea observables (that are 
integrated over Bjorken x) in comparison with the data. It was also shown [TO] that the rf 
meson proposed in [12] gives an almost negligible contribution to the spin fractions of the 
nucleon at the present level of experimental errors not only because of its large mass but 
also due to the small singlet chiral coupling constant. It is therefore ignored in the following, 
and this is consistent with the understanding that, due to the axial anomaly, the rj' meson 
is not a genuine Goldstone boson. 

A major objective here is to construct proton structure functions so that, when integrated 
over Bjorken x, they reproduce the recent successful chiral field theory results for the spin and 
flavor fractions of the proton []10 ||. This goal can be achieved, but only upon neglecting 
helicity non-flip chiral transitions in the polarized splitting functions. The latter 
are also neglected in the unpolarized splitting functions, but this approximation is not nearly 
as critical. 

Standard quark models like the nonrelativistic NQM or light-cone quark models are at 
scale Aqcd, and it is not clear how these models should be evolved to the scale A x in 



the non-perturbative regime where initial quark distributions are needed. In [T3| ratios of 
structure functions were found to be fairly insensitive to the uncertainties inherent in the 
evolution and were therefore recommended as a tool for testing quark models for missing 
physics. An example is the connection between the negative slope of the ratio of unpolarized 
structure functions, -F^/Ff, and a spin dependent attraction between the up and down 
quarks in the proton. In order to construct initial valence distributions, we turn to a more 
successful alternative which has recently been developed in (JT3J, called BBS henceforth. 
This approach is based on quark counting rules or, more precisely, it implements the leading 
power-law behavior at x — > 1 of the helicity-dependent distributions that is known from the 
minimally connected perturbative diagrams in conjunction with constraints from the SU(6) 
proton spin-flavor wave function of the NQM. The initial parton distributions provided by 
BBS have been found surprisingly successful when evolved to the scale of the deep inelastic 
scattering (DIS) data from about 2 GeV |I5] . Therefore, this approach is adopted and used 
as a benchmark to assess chiral dynamics for its ability to generate intrinsic quark- ant iquark 
excitations from initial valence quark distributions of the proton. We shall see that chiral 
field theory then predicts reasonable sea quark distributions. 

The chiral field theory of parton distributions is described in Section II. The antiquark 
and polarized quark distributions are presented and discussed in Section III. Moments are 
given in Section IV and used in Section V to impose quark counting rules from pQCD on 
the valence part of the initial parton distributions. 



II. QUARK DISTRIBUTIONS FROM CHIRAL FIELD THEORY 

If the spontaneous chiral symmetry breakdown in the infrared regime of QCD is governed 
by chiral SU(3)l X SU(3)r transformations, then the effective interaction between the octet 
of Goldstone boson fields $j involves the axial- vector coupling 
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£>int = -^f E^7 M 7 5 A^g (1) 

^Jtt i= i 

which can flip the polarization of quarks: g^ — > g-f + GB, etc. In Eq. [I], the Aj, (z = 
1,2,..., 8) are Gell-Mann's SU(3) flavor matrices, and qa is the dimensionless axial vector 
quark coupling constant that is taken to be 1 here, while 

gf = Au - Ad = A 3 = T + V = (G A /G v ) n ^ p , (2) 

is the isotriplet axial vector coupling constant of the weak decay of the neutron, and Am, Ad 
and As stand for the fraction of proton spin carried by the u, d and s quarks, respectively. 
They are defined by the matrix elements of the singlet, triplet, octet axial vector currents, 
Aff for i=0,3,8 of the nucleon state at zero momentum transfer. It is also common to define 
the hypercharge or octet spin fraction A 8 and the total proton spin 2S Z = AS in the infinite 
momentum frame as 

A 8 = Au + Ad-2As = A£ = Aw + Ad + As. (3) 

The spin fractions that we discussed so far are related to the first moments of the polarized 
structure functions g±(x), g™(x) of the proton which have been measured in a series of SLAC 
(E142, E143, E154 §) and CERN (EMC, NMC, SMC 0) deep inelastic lepton scattering 
experiments with increasingly sophisticated analyses based on pQCD. 

Returning to chiral field theory that we use for a description of DIS structure functions of 
the proton, note that, despite the nonperturbative nature of the chiral symmetry breakdown, 
the interaction between quarks and Goldstone bosons is small enough for a perturbative 
expansion to apply. 

Writing only the flavor dependence of the chiral interaction the SU(3) symmetric Eq. [I] 
has the matrix form 

Lint = §j- EgA^g, (4) 

*Jir i= i 



1 8 

— 1= y Ai$i 



( 4^0 + n + K + \ 



7T~ 



V 



K~ 



-fa ) 



(5) 



From Eq. [I] the following probability densities f(u^ — > tt + + d^),... for chiral fluctuations 
of quarks can be organized as coefficients in chiral reactions: 



d 



-> fd^n-u(n + «J + fd^d^(V + d l) + fdr^fid^° + d i, 



+ fd^K°s{K° + S{), 



(6) 
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and corresponding ones for the other initial quark helicity. After integrating over transverse 
momentum in the infinite momentum frame, the coefficients in Eq. |] become the polarized 
(minus sign) and unpolarized (plus sign) chiral splitting functions, 

PqbM = I d 2 \LU±,{x, k ± ) (7) 



GB/qK-^J ~ J ^Jq-*GBq' 

that determine the probability for finding a Goldstone boson of mass ijigb carrying the 
longitudinal momentum fraction xqb of the parent quark's momentum and a recoil quark 
q' with momentum fraction 1 — xqb for each fluctuation of Eq. [| 

In (non-renormalizable) chiral field theory with cut-off A x of ref. |J, the unpolarized 
chiral splitting functions are 

p + / v _ 9 a XGB(m g + m q ,) 2 /-W. \(m q - m q ,f - t] 

- j2 ^ J_ Al dt ( t -mi B y ' (8) 

where t = — k\ + xgb^ 2 ' — (1 — x^jm^] /{1 — xgb) is the square of the Goldsone boson 
four-momentum. Next we obtain the polarized splitting functions from the unpolarized 
ones using that the latter contain the sum of helicity non-flip and flip probabilities, while 
the former contain the difference of these chiral probabilities. Since quarks are on their 
mass shell in the light-front dynamics used here, the axial vector quark- Goldstone boson 
interaction is equivalent to the simpler 75 interaction. Except for a common overall factor, 
the relevant unpolarized chiral transition probability is proportional to 

- -*r[(7 • p + 77^)75(7 • p' + 777,3)75] = 2(pp' - m q m' q ) = (m q - m' q ) 2 - k 2 , (9) 



1 1 g 

1 



where 2pp' = m' q + m 2 Q — k 2 . This is the numerator in Eq. |8] which can also be written as 

■[k 2 ± + [m q/ - (1 - x GB )m q ] 2 ] (10) 
J- — x GB 

and has the following physical interpretation. Recall that the axial-vector quark- Goldstone 
boson coupling 7^75/^ in Eq. |l] involves the spin raising and lowering operators o\ ± ia-i 
in a scalar product with the transverse momentum of the recoil quark, which suggests 
that the k^ term in Eq. [X0| represents the helicity-flip probability of the chiral fluctuation, 
while the longitudinal and time components 7s7± induce the non-spinflip probability which 
depends on the quark masses. This can be seen from the helicity non-flip probability 

I -/ 1 2 1 — / |2 / / / \2 / 1 /1 1 \ 

K^5 m tI = K75«|| ~ (m q -xm q ) , x = 1-xqb, i 11 ) 

using light cone spinors and suppressing the spinor normalizations. Thus Eq. |TT] identifies 
the mass term in Eq. |10| as the helicity non-flip chiral transition. Similarly, the helicity-flip 
probability is obtained from 

K75«tI 2 ~ Pi + ^' 2 p1 - z'Pi ■ Pi (12) 

which, in frames where = 0, reduces to kj_, i.e. the net helicity flip probability generated 
by the chiral splitting process. In the nonrelativistic limit, where |p'jj << |pi| << m q , 
clearly non-spinflip dominates over spinflip, while spinflip dominates at high momentum. 



5 



The polarized splitting function P therefore has the same quark mass dependence as 
P + , but involves the helicity flip probability with the opposite sign, i.e. has 



1 - x G b 

in its numerator, which agrees with ref. fl6[ . 



h k ± + l m q' ~ (! - X GB )m 

Xnn 



The q — > q'GB splitting functions are expected to satisfy number and momentum sum 
rules . For example, if the quark q is accompanied by two spectator quarks of a nucleon 



which splits (virtually) into a Goldstone boson and an octet baryon B containing the recoil 
quark q' along with the spectator quark pair, then the splitting function in Eq. [13] becomes 
the hadronic f GB /N(x GB ) for N —> B + GB, which equals f B /N^ — x GB ), since the recoil 
quark q' has longitudinal momentum fraction 1 — x GB . Therefore, there are equal numbers 
of Goldstone bosons emitted by the nucleon, (ti)gb/n, an d virtual baryons accompanying 
them, (ti)b/n- This symmetry is violated by Eq. |8] of the model 0. Furthermore, momentum 
conservation holds, viz. (xgb)gb/n + (xb)b/n = ( n ) G B/N- In order to satisfy the symmetry 
and momentum conservation laws, we follow ref. [|T7| and replace the cut-off in Eq. ^| by a 
vertex function G in the splitting function integrand 

f ± f ,2 \ _ 9GBgg> \G G B qq '{x GB ^\)\ 2 ((1 - x GB )m q - m q i) 2 ± k| 
W(^,kJ- lg7r3 Xgb{1 _ Xgb) {l- XGB ){ml-M GBqi y ' {U) 

where is the transverse momentum of the recoil quark q' whose longitudinal momentum 
fraction of the parent quark's is 1 — x GB . Note that 

f ± = f ± = f ± = f ± 

The vertex function G G B qq ' in Eq. takes into account the extended structure of the 
Goldstone boson-quark system. It is taken here to depend on the inverse of the same 
covariant propagator {m 2 q — M GBql )~ l that occurs in the splitting function in Eq. [TB|, 

n ( 12X i m \ ~ M GB q >(x G B,kj) \ 

G GBqq >(x GB , kl) = exp I — ^ I (14) 

with the invariant mass squared of the quark-Goldstone boson system GB + q' 

MS B ,W|) = I?k±2i + ?M. (15) 

x G b i — ^GB 

Such models satisfy the symmetry constraint mentioned above. 

The parameter A in Eq. [14] controls the size of the Goldstone boson-quark system so 
that 2/ A ~< t 2 > GB qi ^ ^ fm. If we relate the monopole shape of niq^p 2 ^ mentioned in the 
Introduction to a Gaussian shape, then A ~ 2^/5Aq GD /0.78 ~1.14 GeV is obtained. Such 
hadronic descriptions of meson clouds of baryons are currently studied by many groups. 
Since the vertex-times-propagator form of the splitting functions in Eq. [13] admits a wave 
function interpretation for the virtual hadronic process — > B + GB, a light-front wave- 



function parametrization has been adopted in ref. [18| for — > A + K. 
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The pion-quark coupling constant in Eq. [13] may be obtained from the Goldberger- 
Treiman relation (GTR) 

(3) m iV na \ 

QttNN = 9 a ~~n ' ( 16 ) 

fir 

where f w ~ 93 MeV is the pion decay constant, in conjunction with the NQM coupling 
constant relation 

9-nqq 1 = -9-kNN- (17) 

If a constituent quark mass instead of the nucleon mass is used in a GTR for constituent 
quarks as in Eq. ||, 

gnqq' = 9 A —, , (lo) 

which is common in chiral quark models at scale Aqcd, then the factor 3/5 in Eq. |D] 
is replaced by qa/^9a '> where qa is the quark axial- vector coupling constant often taken 
to be 3/4 or 1. Since the factor 3/5 in Eq. [L7| is scale (and model) dependent, we set 
g-irqq' — a x9A ^ m N I fix and vary the chiral strength parameter a x within the limits | < a x < |. 
In the SU(3) symmetric case fx ~ fn for the kaons, and we set gxqq' = g^qq 1 - 
The integral over Bjorken x of the splitting function 

PcB/q = \a G B/Q\ 2 = J o dx J d 2 k ± f^ GBql , (19) 

is the probability of finding a virtual Goldstone boson GB in the initial quark q. This in 
turn is related to the Fockspace expansion of the dressed quarks (in light front field theory 
where it is well defined), e.g. 

\U) = Vz\uo) + a n/u \dir + ) + -a^/ulwi ) + -a v/u \ur]) + a K/u \sK + ) + (20) 

where the factors p n o = 1/2 for the neutral pion, p v = 1/6 and = 1 = px± = 1 come from 
the flavor dependence in Eq. |5|, Z = 1 — P u is the wave function renormalization constant 
and all masses and coupling constants are renormalized quantities. 

From the u and d quark lines in Eq. ^ the total meson emission probability P u of the u 
quark is given to first order in the Goldstone fluctuations by 

Pu = f Q dx J d 2 k ± [i/+_w> + lfLu V + fL d , + + fLsK + ] =T,P m f Q dxP^ qm , (21) 



q,m 



while the total strange quark probability density 



R = / dx I d 2 k 



2 



(22) 



can be read off the s quark line in Eq. ^| and total meson emission from the d quark is equal 
to that from the u quark, Pa = P u - 
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The proton's probability composition law P,i~0|,ll,|6f 

(1 - P M )(wJw T + u\u i + dfd T + djd;) 



+ 



(23) 



q,m 



q,m 



q,m 



q,m 



is a useful tool in conjunction with Eq. || for keeping track of chiral fluctuations from valence 
quark distributions q°(x). The latter will be addressed in Sect.V. 

The splitting of quarks into a Goldstone boson and a recoil quark corresponds to a factor- 
ization of DIS structure functions that is equivalent to a convolution of quark distributions 
with splitting functions at a scale to be discussed in Sect. III. From Eqs. |6|j23| the unpolar- 
ized quark distributions off valence quark distributions, q°(x), and up to second order chiral 
fluctuations, are given by convolution integrals of the general form 

1,1 dx\ 



Jx JL\ JL\ 



q ®p + ® q °= ! ld ^<iGB{-) r —p^GsC-m^) 

Jx X\ X\ Jx\ X<i Xi 



x , r 1 dx<} 



so that 



u(x) = Zu° + P 



d—fun 



® If 



+ <+ ® PLsK + ® «° + ® ® («° + d°) + ® P 9 t w ® («° + d°) 

+ l<^®P q tq 7T0 , v ®(u°-d ), (24) 
d(x) = Zd° + P+^ + ® u° + Qp+^o + -P^) ® d° + ® ^ ® 
+ <o ® P+^o ® + \€o ® P 9 t^ ® («° + d°) + ® P q t qv ® (u° + d°) 

-\<^®P^, v ®(^-d% (25) 

* = P u^sK + + P+^ ® d° + is° ri ® P q t qv ® (u° + d°). (26) 
Here the q ( Q B are the valence quark distributions of the respective Goldstone boson. 
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Moreover, the u (or d) amplitudes from the 7r° and rj meson exchanges add coherently to 
the total antiquark probability amplitude of the proton. We know from ref. 151,101,111 that 
these interference amplitudes, i.e. the last terms in Eq. |24],^5], cannot be neglected. Here we 
have approximated the (tt°, 77) interference terms by using equal ir° and 77 splitting functions 
and quark distributions q^o « q v « q in these mesons so that 

/^^(^(^/^^(^^(^i^)^ (27) 

~L^)jju^yQ^ (28) 

takes the form of a conventional parton distribution. 



III. ANTIQUARK AND SPIN DISTRIBUTIONS 

In chiral dynamics, antiquarks originate only from the Goldstone bosons via their quark- 
antiquark composition, viz. 

1 1 

|tt°) = -j=(uu - dd), \tj) = ^(uu + dd- 2ss), \K + ) = us, (29) 

etc. Therefore, antiquarks are unpolarized. Small antiquark polarizations are consistent 
with the most recent SMC data M so that we expect only small corrections if we use = u± 
in the spin fractions Am = — + — ttj_, etc., i.e. As = As sea , Au = Ad = As = 0. 

An antiquark in a Goldstone boson which is struck by the virtual photon in DIS is linked 
to a valence quark of the proton by Goldstone boson exchange. Therefore, antiquarks appear 



only in second order of chiral fluctuations. Thus, from Eqs. |q]23| the antiquark fractions of 



the proton are given by double convolution integrals of the general form q ® P q ^, q GB ® q°, 
u(x) = ul- ® P+__ ® d° + ® P^ qn0 ® (u° + d°) + ±u° v ® ® (1* + d°) 

+ lu° 7r0tV ®P q t q ^ 71 ®(u°-d°), (30) 

d( x ) = dl + ® P^ d7T+ ®u° + ® P+_^c ® (u° + d°) + ® P+_^ ® («° + d°) 

-^® P ™o i?? ®(« -rf ), (31) 

s - = ® P+^ ® d° + g0, + ® p+^ + ® M o + Igo g, g, {u o + d o } (32) 

Again, we have approximated the (tt°,i]) interference terms, the last expression in 
Eqs. pO| , pT| , by using equal tc° and r\ splitting functions and the same antiquark distribu- 
tions in these mesons. 
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The s(x) and s(x) distributions of Eqs. |26], [32], respectively, are no longer equal unless the 
s valence antiquark distributions in the K mesons, s° K0 K+ ~ 6(1 — x), are approximated by 
the elastic limit. When spectator quarks are included in addition, but second order chiral 
fluctuations are neglected, then one recovers from Eq. ^ the strange quark distribution, 



s(x), of some hadronic models [17 , 18 



Upon generalizing the formalism of JlO| also to the polarized quark distributions, the 
probabilities displayed in Eq. 15] and Eq. E3] yield 



qi ( X ) = (l-P q ) q °( X ) + J2PmPl 



*qm 



(33) 



m,q 



where the constants p m are defined below Eq. ^ and m denotes a Goldstone boson. The 
corresponding result holds for the other quark helicity. The ellipses in Eq. denote double 
convolution terms with q® from a Goldstone boson m that cancel in — The opposite 
quark helicity on the rhs of Eq. |33] means that helicity non-flip chiral transitions are neglected 
and implies the negative sign of all chiral contributions to the spin distributions 



Aq(x) = (1 - P q ) Aq\x) - £ P mP^ qm ® Aq'°. 



(34) 



m,q' 



The general result involves the polarized splitting functions 



Pq^q'GBiv) = J d2k ±fq^q'GB(y^ k J 



and has the same form as Eq. [34], 

-i dy 



Au(x) = (1 - P u )Au\x) +J x J [Au°(y) (-P„o + -P^ uv ) + Ad°(y)P d 



(35) 



(36) 



Ad(x) = (1 - P d )Ad°(x) + f 1 ^ 

Jx y 



^{y)(\p d ^ + l -P d ^ dv 



(37) 



As(x) = / — 



idy 

x y 



(38) 



except for the replacement of —P + by the corresponding polarized splitting function P~ . 

Due to the minus sign in the term of the polarized splitting functions in Eq. [L3| 
which corresponds to the subtracted terms in Eq.|34], the general reduction of Ag° by chiral 
fluctuations in lowest order is the crucial property responsible for the success of chiral field 
theory for the proton spin fractions. A comparison of Eq. [£| with Eqs.|3B],|37J, shows that 
only the unpolarized splitting function occurs in the polarized and unpolarized spin-flavor 
fractions of refs. P,[H]|-4T2|,|6[ , and the negative sign of P ~ — P + is critical for the reductions 
of the spin fractions by the chiral fluctuations and the remarkable success of these models. 
Thus, when these Aq(x) are integrated over Bjorken x, the lowest moments (which factorize 
into products of moments as we discuss in the next section) reproduce precisely the structure 
of the results for the spin fractions [|H|P] 
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Au = (1 - P)Au° 
Ad={l- P)Ad° 
As 



5 A " 0+Ad °n i+ ^); 



-(Ha/ + a»> i - * 



(39) 



The pure spin flip approximation is valid at high quark momentum or, upon comparing 
P~ = —P + of Eq. |34] with the numerator of the splitting functions in Eq. [13], when the 
quark masses are negligible compared to the transverse quark momentum < |k|j_ >, i.e. 
current quarks are the relevant degrees of freedom. The factor 2/3 in Eq. ^9| comes from 
adding | + | in Eqs. §^,0 using the same splitting function for the 7r° and r\ mesons because 
the SU(3) breaking along the Ag direction (which is the same as for the standard baryon 
and meson mass splittings) does not suppress the rj contribution. In particular, taking the 
elastic limit for the quark distributions in the mesons, q(x 2 ) ~ 6(1 well as the 

elastic-limit delta function for the valence quark distributions q°(xi), but keeping all second 
order chiral fluctuations, yields all our previous integrated spin-flavor fractions [0,^] and 
Eq. ^9] in particular. Comparing the unpolarized quark and antiquark fractions already 
yields the identification 



a(l + 



dxP„ 



all - 



3L) 

V3 J 



dxP: 



sK+i 



(40) 



where a is the chiral coupling strength and e the SU(3) breaking parameter of [10, B|. For 
other chiral models such as |TT| , p!2| similar identifications can be made. 

This detailed comparison indicates that the success of these chiral models can be at- 
tributed to parametrizations that rely on pure helicity flip chiral transitions. Since the 
helicity non-flip chiral transition probabilities depend on the quark masses, these fluctua- 
tions are negligible only when the quark masses are small compared to the relevant transverse 
momentum range. This is not the case for constituent quarks. Therefore, chiral quark mod- 
els for the spin fractions have to include helicity non-flip fluctuations. When this is done as 
in ref . 



16] then the spin fractions disagree with the data. 



IV. MOMENTS 

The proton spin observables Ag are the lowest moments of, or integrals over, the po- 
larized structure functions. Also, the constraint imposed on the quark distributions by the 
momentum involves the first moment Jq 1 xq(x)dx of the quark distributions. Therefore, in 
the context of DIS and chiral dynamics, it is mandatory to consider moments of convolutions 
of quarks or antiquarks with splitting functions. 

It is straightforward to verify that the integral of a convolution, e. g. 

/ q®P®qf(x)dx= \ q(x)f(x)dx- \ P(x)f(x)dx- \ q(x)f(x)dx, (41) 
Jo Jo Jo Jo 



11 



factorizes provided / is a multiplicative function, i.e. f(xy) = f(x)f(y), such as the power 
law x n that defines the n'th moment. In other words, the moment of a convolution is the 
product of the moments of the components of the convolution. 

Applying this general observation to the quark, antiquark and spin distributions of 
Eqs. p4| , p6| . |30| . p2| . |36| . f??].^ yields for the moments q n = Jq x n q(x)dx the following relations 



u n (x n u -\- b n d , d n ci n d -\- b n 



u. 



, _ p+,n , p+,n JO , Z„0 p+,n ( , M \ 



(42) 



1 



1 



Un = U°^ n P^-dn + - A <»,nPTK + ft) + - a ul n P+' n {ul + ft) + -«°P + '»(tl° - ft] 



36 



"n 1 



1 . 

6^' 



(43) 



where the P + > n are the moments of the unpolarized splitting functions and the q®, s° are 
those of the proton and K meson valence quark distributions. The coefficients a n , b n are 
given by 



/ - 7 4- - P + ' n 4- - P+' n 4- ii° P + ' n 

2 7r ° 6 n ,n^u-*dn+ 

-l-I,, P+' n 4- J-7,° P+> n 4- Ii, P+> n 



36 



6 



(44) 



Z, _ p+,n , 1 p+,n , 1 p+,n _ 1 p+,n 



36 



6 



Similarly, from the Ag in Eq. the spin dependent moments are obtained as 



Aw r 



'1 - P n )AvS. 



Ari° Ad! 



/° 



(45) 



(46) 



Ad n = (1 - P d )Ad° 



1 1 

P~' n 4_ p->« 

2 *° 6 77 



A< + P.2» A«° 



(47) 



A*n = ^jc + A< + P d T^oA< 



(48) 



to be compared with Eq. [39| for n = 0. 

As a consequence of the factorization of moments of convolutions, the equality of the 
integrated s = s of the proton follows immediately from Eqs. 42,43 using Sq = 1 for the 



valence s antiquark fraction in the K mesons. In order to avoid more parameters and 
complications we shall use the elastic delta function approximation for the valence quarks 
and antiquarks in the Goldstone bosons in the numerical evaluations. 
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V. QCD CONSTRAINTS AND RESULTS 



There are numerous attempts to describe DIS structure functions based on quark mod- 
els. These quark models provide initial quark distributions at the scale Aqcd as a rule. For 
the NQM [|l|J and a light-cone quark model (LCQM |13| ) the up quark distributions are 
exhibited in Fig.l. The subsequent Q 2 -evolution of structure functions follows from pertur- 
bative QCD, but relies on such initial shapes which contain the non-perturbative aspects of 
the bound state. The low-energy scale of a quark model is often treated as an adjustable 
parameter which is usually too small for a perturbative evolution to be valid. It is fair to 
say that, even with chiral fluctuations, quark models fall short of explaining the deep inelas- 
tic structure functions and the spin of the proton [16[], quite apart from rather significant 
uncertainties in their evolution from Aqcd to about A x . 

The success of chiral fraction models based on pure helicity flip chiral fluctuations sug- 
gests that we construct initial quark distributions for the case where current quarks are 
the relevant degrees of freedom. Upon comparing with the momentum dependence of the 
dynamical quark mass discussed in the Introduction, the relevant momentum scale is about 
| GeV to A x rather than Aqcd of the NQM. It is important to realize in this context that 
the use of the spin fractions Ag° from the spin-flavor proton wave function of the NQM 
in Eq. RDI is not only justified at scale Aq C d but at A x as well because the (for As = 0) 

(3) (S) 

equivalent isotriplet and octet axial charges g\ ,g\ are conserved in the chiral limit and 
scale independent. 

Instead of the NQM we use the BBS approach to construct valence quark distributions. 
In contrast to BBS, we impose the quark counting rules directly on the valence quark dis- 
tributions. There are no parametrizations for strange quarks or antiquarks which we expect 
to be generated by chiral fluctuations. It is an immediate consequence that initial parton 
distributions where the strange quarks contribute nothing (for some examples, see |2t|) are 
ruled out by chiral field theory. 

The constraints 



1 rl 

u°(x)dx = 2, d° = / d°(x)dx = 1 (49) 



JO 



on valence quark numbers are characteristic features of quark models, such as the nonrela- 
tivistic NQM and light cone quark models. Moreover, the quark distributions satisfy q° — > 
as x — > in these models. Hence for illustration purposes only, we parametrize such valence 
quark distributions q°(x) in terms of the shapes x a (l — x) 13 with a common positive expo- 
nent a and f3 =3, or 5. As a consequence of the positive exponent, fewer non-leading terms 
are needed, which saves more than half of BBS's parameters, viz. 

u^{x) = Ax a {l - xf, = Cx a {l - x) 5 , 



d^(x) = -x a (l-x) 3 , dj(x) = C'x a (l - a;) 5 , (50) 

where the anti-parallel helicity quark distributions are suppressed by the factor (1 — x) 2 
and the exponent (3 =3 for it? and d® reflects the quark counting rule for the proton that 
guarantees the correct falloff of the leading form factors with q 2 . The choice A/ 5 for the 



13 



leading coefficient of c?° in Eq. |5y is motivated by the SU(6) structure of the NQM and its 
valence quark distributions for large x. When the other ratios of coefficients are also chosen 
as in the SU(6) spin- flavor wave function of the proton (cp. Eq.13 of 
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A 9 
C=t C' = Z -A, 
5 5 



then we obtain 



A = gfm =l (j + a)/6 d 



a 



161 36Ad° 



1/2 



11 

y 



i 

2 



(51) 



(52) 



for an initial value Ad° ~ —0.2 which increases to -0.4 by chiral fluctuations. The corre- 
sponding up quark distribution with chiral fluctuations is close to the dot-dashed line in 
Fig.l, where the exponent a = 1/2 and Au° = 4/3 are chosen to simulate the NQM. It is 
straightforward to express the coefficients C, C in Eq. [SO for that case in terms of A and 
a from the two constraints of Eq. The formulas for the coefficients in Eq. 50 are given 



in the Appendix and their numerical values are displayed in Table 2. Upon imposing the 
observed values for , \ or equivalently Au°, Ad° as well, the relation Am — 5 Ac? = 3 
is obtained, which is satisfied by the NQM values Am = 4/3, Ac? = —1/3. Thus, for fixed 
values of the exponent a and Am , which are not adjusted to any observable, the chiral 

(3) 

strength is adjusted to fit g A with chiral fluctuations, yielding the results in Figs. 1,2. The 
case with a smaller Au° = 1.1 is supposed to simulate a relativistic quark model. The price 
one pays for selecting a positive x-power is that the resulting quark distributions shown as 
dot-dashed line in Fig.l for exponent a = 1/2 and solid lines in Fig. 2 become unrealistically 
small for x — > 0, which is a typical quark model feature. Clearly, these quark distributions 
have a lower peak than those of the NQM and LCQM, and it is shifted to lower x values, 
features that make them look like evolved from quark models in Fig.l at Aqcd to a higher 
scale. The comparison of the dotted curve in Fig. 2 which is the associated valence u°(x) 
quark distribution without chiral fluctuations with the upper solid line shows the typical 
size of chiral fluctuations; they shift the peak to lower x but increase its magnitude for the 
u quark, while it is lowered for the d quark. 

The integrated spin observables obtained with chiral fluctuations are given in Table 1 
for the cases shown in Figs. 1,2, 3. Note that for the quark model simulations with exponent 
a =1/2, despite fitting only one parameter, the chiral strength a x , to the axial charge of 
the nucleon, most of the other spin observables are not too far off. In Fig.3, though, the 
exponent a =0.1 is deliberately chosen small to get closer to the BBS result for u(x). For the 
smaller exponent a = 0.1, the total quark momentum constraint u\ + c?° = 0.52 is imposed 



instead of the quark number constraints of Eq. W, which also helps getting (the dotted line 
in Fig.3) closer to the BBS (dot-dashed) result, except at small x. For this case, nearly all 
integrated spin fractions in Table 1 are already in good agreement with the data, a success 
for which chiral fluctuations are essential in conjunction with current quark masses. Clearly, 
all these quark distributions become unrealistically small for x < 0.1. 

When the light up and down quarks are taken to be constituent with m q = m^/'i = 313 
MeV, and the strange quark with m s = mj — mjv + m q = 567 MeV, then most numerical 
results become worse; the chiral strength has to be increased, while the proton spin AE 
moves up to about 0.5, confirming the results of ref. [16 
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Since the recent HERA data |2T| clearly show the rise of the unpolarized F%(x) at small 



x, which the quark model-type distributions in Figs. 1,2, 3 fail to reproduce, in the next step 
we incorporate Regge behavior imposing a negative common (isoscalar) power x~ aR on the 
valence distributions. Just like BBS, we find that more non-leading terms must be introduced 
to make the proton spin fractions Aq and the Gottfried sum rule (GSR) convergent. The 
minimal Regge-type parametrization is given by 

v°(x) = x~ aR [A + B{l -x)](l -x) 3 , uj(x) =x~ aR [C + D{l -x)](l - x) 5 , 



(1-x) 3 , d\{x) = x- aR [C + D'{l-x)]{l-xf. (53) 

For Aw°(x), Ad°(x) to be finite for x — > 0, the constraints 

A + B = C + D, - + B' = C + D' (54) 
5 

are required. For the Gottfried sum rule to be convergent, or isospin symmetry for x — > 0, 
requires 

A + B + C + D = 4 + B' + C + D'. (55) 
o 

Moreover, we impose fixed values for Au°, Ad° and the GSR, or equivalently u® — dg, and 
the total quark momentum u\ + d\. Again, it is straightforward but tedious to determine 
the coefficients in Eq. |53| from these conditions; they are given analytically in the Appendix 
and numerically in Table 2. 

Chiral fluctuations will not change the Regge exponent in the dressed quark distributions 
because of the asymptotic behavior of convolutions, e.g. from q°(x) ~ x~ aR we obtain 

P ± ®q° ~ x~ aR f 1 dtP^t- 1 -"*, x -> 0, (56) 
Jo 

where the moments of the splitting functions remain finite because the Gaussian in Eq. |14] 
goes to zero faster than any power at the endpoints. 

We adopt the Regge exponent a# = 1.12 of BBS. Then we adjust the chiral strength 
to the spin AS of the proton roughly, vary Au° — Ad°, Au° + Ad° (which changes AS 
a bit) to match the observed values for g\, g\ with chiral fluctuations, impose u® — g?q to 
match the GSR and v![ + d\ w 1/2 of the total quark momentum. The results are displayed 
as solid lines in Figs. 3-8. The d quark distribution in Fig. 4 deviates more from BBS, but 
the s quark distribution in Fig. 5 (solid line) predicted by chiral fluctuations is again close 
to the BBS (dot-dashed) curve, whereas the quark model case with exponent a =0.1 is 
much too small at small x. The spin dependent Ag distributions in Figs. 6-8 show a similar 
pattern. The polarized proton structure function shown in Fig. 9 is crudely evolved by a 
multiplicative 1 — a s /7t ~ 0.92 to the data at — q 2 =10 GeV 2 as in fll4| . The ratio of 
unpolarized F^jF^ — > 3/7 as x — > 1 with and without chiral fluctuations. Finally, the 
antiquark ratio u/d in Fig. 10 is not even approximately constant. 



x 



-<*r 



A 



+ B'(1 
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VI. CONCLUSIONS 



A comparison of the chiral field theory of quark structure functions with previous chi- 
ral models P,|l0Hl"2],|6f for the spin fractions Aq shows unambiguously that the remarkable 



success of these models is due to parametrizations that correspond to neglecting helicity 
non-flip chiral fluctuations in the (integrated) polarized splitting functions which are sizable 
for constituent quarks. Thus, standard chiral dynamics can explain the spin observables 
of the proton only for the case where the dynamical quark masses are already small. Fur- 
thermore, the present analysis sheds light on the failure of standard chiral dynamics at the 
scale Aqcd involving constituent quarks |TjJ whose larger masses significantly reduce the 



negative chiral contributions to the spin fractions. When quark distributions are constructed 
according to the quark counting rules of pQCD in conjunction with constraints from the 
SU(6) spin-flavor wave function of the proton, the GSR and total quark momentum fraction, 
then chiral dynamics generates sea quark distributions of correct magnitude and shapes and 
predicts reasonable strange quark distributions in particular so that the spin fractions of the 
proton also agree with the data. 

While chiral dynamics with the (ad-hoc) pure spinflip prescription appears to be success- 
ful in understanding the proton spin fractions, its justification in the context of chiral quark 
models may be provided by instanton dynamics of zero- mode quarks. In fact, 't Hooft's 



interaction maintains pure left- to right-handed quark transitions, while light quarks pick 
up a dynamical mass upon propagating in the instanton background field. J^J However, 
axialvector coupling of constituent quarks with Goldstone bosons in standard chiral field 
theory is inconsistent with the proton spin fraction data. 

VII. APPENDIX 

The coefficients of the valence quark distributions for the cases of Fig.l (dot-dashed line) 
and Fig. 2 (upper solid and dotted lines) with a =1/2 are given by 

, 2 + Au° 4 .. . „ 2-Ait° .. 

a = 12 nj =1 (j + a), c= 24Q ng = i(,7 + a), 

8-A^° 

1200 J ~ 1KJ h V ' 

those for Fig. 3 (dotted line) with a =0.1 and total quark momentum constraint are 

A = + d?)(5 + a) + (Am + Ad°)(l + a)]U^ =2 (j + a), 



A A?; 

C = — (6 + a)(5 + a) - —itfJj + a) 
20 V A 1 120 ' 



A AH 
C = 100 (6 + a)(5 + a) " 120 "5= l(j + a) ' (58) 
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while those of the Regge (solid lines in Figs. 3-8) and BBS (dot-dashed lines in Figs. 2-8) are 
given by 



A=- (Au° - Ad° + u° - d°) U 5 j=2 (j - a R ), 



-405 



[(12 - a R ){2 - a R ) + (7 - a R )(8 - a R ) + 30] = [6(6 - a R )(7 - a R )(8 - a R ) 



16(7 - a R )(8 - a R ) + 1200 + (2 - a R )(732 - 118a fl + 6a 2 R )} 
1. 
6 1 



+ «nj =3 (j - a R ) [K + d°)(8 - a fl ) + (Am + Arf°)(2 - a*) 



(59) 



Am 

40C = - 



40C = ^nj =2 (j - + 85(12 - a R ) 



+ ^[4(7 - a R )(Q - a R ) + 732 - U8a R + 6a£], (60) 
o 

Ar/° 

40C" = nj =2 (j - «i?) + 85(12 - 

+ 4[732 - 118a fl + 64 - 4(7 - a R )(Q - a R )\. (61) 
o 



B' = -A + B, D' = A + B-C. (62) 
5 
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Table 1 Quark Spin Fractions of the Proton for Figs. 1-8 
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Table 2 Coefficients of Valence Quark Distributions in Eqs. |50| , [53 
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BBS 


141 




cx^ — u.oy 


0,^ — u.o 




ru fl A 

U>X U.4I: 




Figs.2-8 


n, — 1 


n, — 1 
a - 2 


a — u. 1 


ru„ 1 IO 

Oj? — 1-1^ 








An — 4 
— 3 


A,,0 _i 1 












Fies 1 2 


Fig's 2 


Figs. 3-8 


Figs. 3-8 


A 


3.78 


4 


16.406 


15.258 


8.091 


5.565 


B 


-3.672 








-5.551 


B' 


-0.645 








-1.099 


C 


2.004 


5.866 


7.918 


-0.354 


-1.047 


c 


3.230 


11.730 


12.141 


2.974 


0.920 


D 


-1.892 








1.061 


D' 


-3.118 








-0.906 



19 



FIGURES 




FIG. 1. Quark distribution xu(x) of NQM [18] (solid), LCQM [11] (dashed), both at A QC d, 
case of 3rd column of Table 1 (dot-dashed) with chiral fluctuations at A x . 
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FIG. 2. Quark distribution xu{x) of 3rd column of Table 1 (top solid) with corresponding 
valence xu°(x) (dotted) without chiral fluctuations, case of 4th column of Table 1 (lower solid) 
compared with BBS [12] (dot-dashed). 
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FIG. 3. Quark distribution xu(x) of 5th column of Table 1 (dotted), Regge case of last column 
of Table 1 (solid) compared with BBS [12] (dot-dashed). 
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FIG. 4. Same as Fig.3 for down quark xd(x). 
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FIG. 5. Same as Fig.3 for strange quark xs{x). 
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FIG. 6. Same as Fig.3 for polarized Au(x). 
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FIG. 7. Same as Fig.3 for polarized Ad(x). 
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FIG. 8. Same as Fig.3 for polarized As(x). 
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